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1.  INTRODUCTION 


1.1  Background 

For  the  past  few  years,  the  Air  Force  Rocket  Propulsion  Laboratory 
(AFRPL)  has  sponsored  a series  of  theoretical  and  experimental  programs  on  the 
retrieval  of  plume  flow-field  properties  by  analysis  of  the  infrared  radiative 
and  absorptive  properties  of  plumeSr  This  report  continues  the  third  phase  of 
study  in  these  programs  by  The  Aerospace  Corporation.  The  first  phase  (Ref. 
1)  was  a study  of  the  classic  problem  of  retrieving  radial  profiles  of  gas 
temperature  and  concentration  in  cylindrically-symmetric,  gaseous  plumes  from 
transverse  emission  and  absorption  (E/A)  profiles  obtained  in  a fixed  spectral 

bandpass.  The  E/A  profiles  are  defined  in  terms  of  the  radial  profiles  of 

/ 

pressure,  temperature,  and  concentration  (pTc  profiles)  by  integral  equations 
of  radiative  transfer.  Retrieval  of  the  jffc  profiles  from  the  E/A  profiles 
involves  a numerical  inversion  of  these  Integral  equations. 

In  this  first  phase  study,  an  inversion  procedure  was  developed  and 
incorporated  into  the  computer  code  EMABIC.  The  inversion  algorithm  is  an 
iterative  Abel  inversion.  The  well-known  Abel  inversion  procedure  is  valid 
for  optically  thin  sources;  for  the  general  case  of  optical  thickness,  an 
iterative  procedure  is  required.  The  code  has  been  applied  to  several  syn- 
thetic and  experimental  data  and  performs  satisfactorily  as  a diagnostic  for 
most  gas-only  plume  problems.  Some  problems  occur  when  the  temperature 


1.  S.  J.  Young,  Inversion  of  Plume  Radiance  and  Absorptance  Data  for 
Temperature  and  Concentration,  AFRPL-TR-78-60,  U.  S.  Air  Force  Rocket 
Propulsion  Laboratory,  Edwards  Air  Force  Base,  California,  29  September 


profile  has  a deep  minimum  on  the  plume  axis  or  when  the  input  E/A  profiles 
are  particularly  noisy,  even  if  they  are  adequately  smoothed.  This  is  an 
inherent  feature  of  Inversion,  however,  and  is  not  restricted  to  the  method  of 
inversion.  A similar  inversion  code  has  been  developed  at  the  Arnold 
Engineering  Development  Center  (AEDC)  (Ref.  2).  Recently,  a random  error 
propagation  routine  was  added  to  EMABIC  so  that  retrieval  error  could  be 
estimated  automatically  from  E/A  measurement  error  (Ref.  3). 

The  second  phase  of  study  was  the  consideration  of  multispectral  inver- 
sion and  the  effects  of  particle  loading  in  tactical  motor  plumes.  In  multi- 
spectral inversion,  retrieval  is  made  on  the  basis  of  how  E/A  spectra  vary  in 
wavelength  for  a fixed  measurement  line  of  sight.  It  was  found  that  this 
inversion  scheme  is  not  applicable  in  the  infrared  on  either  a monochromatic 
or  wide  band  spectral  scale  near  the  exit  plane  for  small  plumes  with  mild 
temperature  gradients,  such  as  those  characteristic  of  tactical  rocket 
motors.  Even  under  ideal  circumstances,  temperature  and  concentration  re- 
trieval errors  up  to  30 % were  encountered.  The  failure  of  the  method  is 
caused  by  the  lack  of  spatial  resolution  inherent  in  the  inversion  weighting 
functions.  Results  of  this  study  are  reported  in  Ref.  4.  Because  this  method 
failed  for  purely  gaseous  plumes,  it  was  never  applied  to  two-phase  plumes. 


2.  C.  C.  Limbaugh,  W.  T.  Bertrand,  E.  L.  Kiech  and  T.  G.  McRae,  Nozzle  Exit 
Plane  Radiation  Diagnostics  Measurements  of  the  Improved  Transtage  Liquid 
Rocket  Injector  Program,  AED  -TR-79-29,  ARO  Inc.,  Arnold  Engineering 
Development  Center,  Arnold  Air  Force  Station,  Tennessee,  March  1980. 

3.  S.  J.  Young,  Random  Error  Propagation  Analysis  in  the  Plume  Diagnostic 

Code  EMABIC,  AFRPL-TR-81-59,  U.  S.  Air  Force  Rocket  Propulsion 

Laboratory,  Edwards  Air  Force  Base,  California,  July  1981. 

4.  S.  J.  Young,  Multicolor  Inversion  Diagnostic  for  Tactical  Motor  Plumes, 
AFRPL-TR-80-30,  U.  S.  Air  Force  Rocket  Propulsion  Laboratory,  Edwards  Air 
Force  Base,  California,  May  1980. 
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It  was  decided  to  revert  to  the  multiposlticii  inversion  diagnostic  of  the 
first  study  phase  and  to  pursue  its  application  to  two-phase,  tactical  rocket 
motor  plumes. 

The  primary  goal  of  the  first  part  of  this  third  phase  was  to  determine 
quantitative  limits  of  particle  loading  in  realistic  tactical  rocket  motor 
plumes  for  which  the  gas  properties  of  the  plume  could  be  retrieved  with 
standard  E/A  inversion  diagnostics  without  having  to  account  for  the  radiat- 
ing, absorbing,  and  scattering  effects  of  the  particles.  The  procedure  was  to 
assume  flow-field  properties  for  two-phase  plumes  of  interest,  generate  E/A 

•ff 

profiles  with  account  of  particles  using  a single-scattering  plume  radiation 
code  ( EAPROF , Ref.  5),  retrieve  the  gas  properties  from  the  E/A  profiles  with 
the  gas-only  inversion  code  EMABIC  under  the  assumption  that  the  profiles  were 
caused  by  gas  alone,  and  compare  the  retrieved  gas  properties  with  the  assumed 
properties.  Generally,  the  degree  of  particle  loading  was  treated  as  a pa- 
rameter. The  work  focused  on  tactical  rocket  motors  where  the  particle  load- 
ing level  is  small,  that  is,  to  motors  where  the  particulate  material  is  added 
to  the  fuel  only  as  a stabiliser  (e.g.,  AI2O3),  or  to  motors  where  the  plume 
particulate  results  from  chemical  reactions  (e.g.,  carbon)  but  not  to  motors 
in  which  the  major  fuel  is  Itself  a metal.  The  procedure  was  applied  to  a 
minimum  smoke  propellant  (MSP)  model,  an  advanced  liquid  propellant  (ALP) 


When  particle  effects  are  accounted  for,  the  A of  E/A  should  more 
properly  be  interpreted  as  "attenuation"  rather  than  just  absorption 
since  signal  reduction  by  particle  scattering  is  included.  In  this 
report,  "extinction"  and  "extinctance"  are  sometimes  used,  along  with 
"attenuation,"  to  denote  the  combined  effects  of  absorption  and 
scattering. 

5.  S.  J.  Young,  User's  Manual  for  the  Plume  Signature  Code  EAPROF,  AFRPL-TR- 
81-08,  U.  S.  Air  Force  Rocket  Propulsion  Laboratory,  Edwards  Air  Force 
Base,  California,  January  1981. 
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model,  and  a reduced  smoke  propellant  (RSP)  model.  The  results  are  reported 
In  Ref.  6. 

Two  Important  results  were  obtained.  The  first  is  that  the  limit  of 
particle  loading  at  which  reasonable  « 10%  error)  retrieval  results  can  be 
obtained  is  generally  smaller  than  the  nominal  loading  level  for  the  plume. 
The  important  implication  of  this  result  is  that,  to  the  extent  that  the 
systems  studied  are  typical,  E/A  diagnostics  on  plumes  generated  by  even  low- 
and  reduced-smoke  type  propellants  require  some  account  of  particle  effects. 
The  second  result  is  that  the  maximum  loading  level  for  acceptable  gas  tem- 
perature retrieval  is  much  higher  (about  an  order  of  magnitude)  than  that  for 
gas  concentration  retrieval.  Consequently,  in  applications  where  temperature 
retrieval  is  of  primary  concern,  the  use  of  gas-only  E/A  diagnostics  may  be 
justified  even  though  the  total  retrieval  results  may  be  substantially  In 
error. 

Analysis  was  also  made  of  a procedure  in  which  first-order  corrections 
were  made  to  the  total  gas-plus-particle  E/A  profiles  by  using  particle-only 
E/A  profiles  obtain* d outside  the  gas  absorption  band.  The  corrected  profiles 
provided  better  estimates  to  the  gas-only  profiles  needed  in  the  gas-only 
inversion.  The  particle  loading  limit  for  valid  use  of  this  procedure  is  the 
value  for  which  the  total  extinction  of  radiation  by  particles  over  a full 
diameter  of  the  plume  is  about  10  percent.  (Note  that  if  this  condition  is 
met,  then  the  condition  is  also  met  that  the  attenuation  by  scattering  alone 
over  this  path  is  less  than  about  10  percent.  The  latter  condition  is 


6.  S.  J.  Young,  Retrieval  of  Flow-Field  Gas  Temperature  and  Concentration  of 
Low-Visibility  Propellant  Rocket  Exhaust  Plumes,  AFRPL-TR-82-13,  U.  S. 
Air  Force  Rocket  Propulsion  Laboratory,  Edwards  Air  Force  Base, 
California,  March  1982. 


required  by  the  single-scattering  assumption  used  in  the  work. ) For  the  two 
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cases  (MSP  and  &SP)  where  the  nominal  loading  limit  roughly  satisfied  this 
condition,  the  use  of  the  procedure  resulted  in  retrieved  gas  properties  that 
were  accurate  to  within  the  convergence  criteria  set  on  the  inversion.  For 
the  ALP  model,  the  nominal  loading  level  was  well  above  this  limit,  and  the 
retrieved  results  were  poor. 

1.2  Scope  of  Present  Study 

The  present  work  expands  the  work  begun  in  Ref.  6 in  two  significant 
aspects.  First,  the  restriction  of  the  first-order,  off-band  (FOOB)  correc- 
tion procedure  that  the  plume  be  optically  thin  to  total  extinction  is  re- 
laxed. This  is  accomplished  essentially  by  using  the  full,  two-phase,  single- 
scattering radiation  model  in  the  inversion  code  EMABIC  as  well  as  In  the 
plume  radiation  code  EAPROF.  The  second  significant  expansion  is  that  re- 
trieval diagnostics  for  particle  properties  are  developed.  The  previous  work 
on  two-phase  plume  diagnostics  assumed  that  all  particle  properties  (e.g., 
species,  index  of  refraction,  spatial  and  size  distribution)  were  known.  The 
current  work  defines  procedures  for  retrieving  the  volume  cross  sections  for 
particle  absorption,  scattering,  and  extinction,  as  well  as  the  particle 
scattering  phase  function  and  radial  temperature  profile.  (Further  analysis 
is  underway,  but  not  reported  here,  on  the  retrieval  of  the  particle  size  and 
radial  distributions  from  the  scattering  phase  function  and  volume  extinction 
cross  section.)  The  plume  emission/attenuation/scattering  (E/A/S)  data  re- 
quired for  application  of  these  procedures  are  discussed  in  Section  1.3. 

As  in  the  preceding  work  of  this  phase,  the  analyses  presented  here  were 
made  with  synthesized  E/A/S  data,  that  is,  data  computed  from  assumed  known 
plume  properties.  No  application  to  experimental  data  was  made.  Also,  since 
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the  main  emphasis  of  this  work  was  on  the  feasibility  of  retrieval  and  not  its 
practical  application,  no  treatment  of  random  or  bias  error  propagation  was 
made.  The  radiative  transfer  theory  used  to  compute  these  synthetic  data  is 
presented  in  Section  3,  and  the  inversion  procedures  for  retrieving  the  plume 
properties  from  the  synthesized  data  are  presented  in  Section  3.  The  entire 
prediction  and  inversion  capabilities  developed  in  these  two  sections  have 
been  written  as  a series  of  computer  codes  and  are  described  in  a separate 
report  (Ref.  7).  Application  of  the  inversion  procedure  was  made  to  the  MSP, 
ALP,  and  RSP  models  described  in  Ref.  5,  Results  and  discussion  of  the  appli- 
cations are  made  in  Section  4,  and  a summary  of  conclusions  drawn  from  the 
results  is  given  in  Section  5. 


1 .3  Emission/Attenuation/Scattering  Data  Requirements 

The  diagnostic  scheme  developed  in  Section  3 requires  traditional  trans- 
verse E/A  data  at  two  wavelengths  and  angular  scattering  data  at  one  of  the 
two  wavelengths  (Fig.  1).  Aj  is  a spectral  position  at  which  both  the  gas  and 

I 

particles  of  interest  radiate  and  absorb.  a position  at  which  only  the 
particles  radiate  and  absorb.  A ^ should  be  as  close  to  A^  possible  so  that  It 
can  be  assumed  that  the  particle  optical  properties  are  the  same  at  the  two 
wavelengths.  The  transverse  E/A  data  required  at  these  wavelengths  are 


N(z) 

t(z) 


N (z) 
P 

rp(z) 


o w.  < R A = Aj  (gas  + particles) 

o < z < R A =»  (particles  only) 


7.  S.  J.  Young,  Us»r ’ s Manual  for  the  Flow-Field  Diagnostics  Code 

EMABIC,  AFKPL-TR-82-37",  U.  S.  Air  Force  Rocket  Propulsion  Laboratory, 
Edwards  Air  Force  Base,  <’a lifornia,  February,  1983. 


where  z is  Che  transverse  coordinate,  R is  the  plume  radius,  N is  the  plume 
radiance,  and  t is  the  plume  transmittance  to  external  radiation.  An  axisym- 
metric,  cylindrical  plume  is  assumed. 

The  scattering  data  are  required  only  at  X In  order  to  get  enough 
scattered  intensity  to  measure  accurately,  the  source  would  probably  have  to 
be  a laser.  This  source  is  scanned  in  the  transverse  direction  just  as  the 
E/A  sources  are  scanned.  Angular  scattering  data  are  required  between  G = o 
and  3 = it  where  the  reference  is  the  z = o line  and  the  plume  center  line. 
The  scattering  function  used  in  this  work  considers  scanning  from  z =*  o to  2 ■= 
R and  angular  measurements  on  the  opposite  side  of  the  plume  (as  indicated  in 
Fig.  1).  The  field  of  view  of  the  scattered  radiation  sensor  should  be  wide 
enough  to  see  the  entire  width  of  the  plume.  Also,  the  distance  from  the 
sensor  to  the  plume  should  be  much  larger  than  R.  Figure  1 indicates  a single 
detector  that  scans  in  8.  A more  probable  arrangement  would  be  an  array  of 
detectors  at  various  0.  In  any  case,  the  required  data  are 


f (z,0) 


0 < 0 < IT 

o < z < R 

X m x2 


where  f is  the  scattered  power  divided  by  the  laser  power.  In  general,  the 
variation  in  9 is  needed  for  each  transverse  position.  A simplified  retrieval 
analysis  is  presented  later  (Section  3.3)  in  which  only  the  variation  at  z * o 
is  required. 
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2.  RADIATION  TRANSPORT  EQUATIONS 


The  radiation  model  for  predicting  E/A  profiles  from  known  plume  proper- 
ties was  derived  and  discussed  in  detail  in  Refs.  5 and  6.  Briefly,  the  model 
treats  the  problem  of  radiative  transfer  in  a coupled,  two-phase  (gas  and 
particles)  medium.  Radiation  transport  for  the  gas  phase  is  treated  with  the 
Malkmus  statistical  band  model  and  the  Curtis-Godson  nonumiformity  approxima- 
tion. Radiation  transport  for  the  particle  phase  is  treated  with  the  single- 
scattering  approximation.  In  this  approximation,  only  radiation  that  is 
emitted  directly  on  and  along  an  observation  line  of  sight,  or  is  scattered 
into  the  line  of  sight  by  a single  scattering  contributes  to  observed  emis- 
sion. In  transmission,  no  radiation  scattered  into  the  line  of  sight  is 
counted  as  transmitted  signal. 

A review  of  this  model  for  a line  of  sight  through  a general  medium  Is 
given  in  Section  2.1.  In  Section  2.2,  the  equations  are  transformed  to  cylin- 
drical coordinates  for  application  to  the  plume  problem,  and  written  in  a form 
suitable  for  both  prediction  and  inversion.  In  Section  2.3,  the  addition  to 
the  model  required  for  calculating  the  laser  scattering  function  f(z,9)  is 
developed. 

2. 1 General  Line-of-Sight  Equations 

For  a line  of  sight  8 through  a general  two-phase  medium  (Fig.  2) , the 
observed  emission  N is 

L 

N - / Ta(s)  T (s)  t(8)  {Qt(s ) + QgCs ) } ds  (1) 

o a 

Qt  (s)  - J (s)  + y(s)  J (s)  (2) 

r O 
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(3) 


J (8)  * a(8)  B (s) 
P P 


J (8)  * <(8)  B (b) 
g g 


Qc(b)  - f Jo  lQ„<S,fl)  + QmCs.fi)]  dn 


S 4tt  S 


Jg(s*fl)  - ^ p(s,n) 


Q (8,n)  - T [L(s,ftJ  N (8,n? 

E a t(s)  e 

L(s ,fi) 

Qt(s ,Sl)  = / T (a)  Q (a, s)  da 

T o a t ( 8 ) T 


Q (a,s)  » J (a)  + y(s  + a)  J (a) 

^ r O 


J (a)  * a (a)  B (a) 
P P 


J (a)  - <(a)  B (a), 
g g 


All  of  the  quantities  appearing  in  these  equations  have  been  defined  in  Ref. 
6,  and  only  a brief  explanation  will  be  given  here.  r^Co)  is  transmit- 
tance between  s * o and  general  position  s effected  by  particle  absorption. 
t (s)  and  t(s)  are  the  corresponding  transmittances  for  particle  scattering 

p 

and  gas  absorption.  Q^,  is  the  thermal  source  function  and  is  composed  of  the 
particle  component  Jp(s)  and  the  gas  component  Jg(s).  a(s)  and  0(s)  are  the 


volume  cross  sections  for  particle  absorption  and  scattering,  respectively, 
and  <(s)  is  the  volume  cross  section  for  gas  absorption.  y(s)  is  a band  model 
function  that  accounts  for  the  deviation  of  band  absorption  from  the  Lambert 
and  Beer  absorption  laws,  and  varies  with  assumed  spectral  lineshape  and 
approximation  used  to  treat  the  optical  path  nonuniformity.  Bp(s)  and  B^(s) 
are  the  Planck  function  evaluated  at  the  particle  and  gas  temperatures  at  s, 
respectively.  Qg  is  the  scattering  source  function  and  is  an  integral  over 
all  directions  of  radiation  arriving  at  s from  the  scattering  line  of  sight 
o.  ft  is  the  scattering  solid  angle  described  by  the  scattering  angle  0 and 
azimuthal  angle  <j>.  The  integration  weight  given  to  each  scattering  line  of 
sight  is  Jg(s,ft)  and  is  directly  proportional  to  the  particle  scattering  phase 
function  p(s,ft).  Qg  is  the  contribution  to  Qg  from  external  radia- 
tion Ng(s,ft)  incident  at  the  boundary  position  o - L(s,ft).  and  t are 

particle  and  gas  transmittances,  respectively,  for  the  indicated  path  lengths. 
QT(s,fi)  is  the  thermal  contribution  to  Qg  from  the  scattering  line  of 
sight.  All  of  the  quantities  comprising  QT  are  as  previously  explained, 
except  that  they  are  written  for  the  line-of-sight  variable  o rather  than  s. 

The  transmission  through  the  medium  along  the  observation  line  of  sight  s 
is  given  by 

L 

-In  t *■  / [y(s)  + k(s)  y( s ) J ds  (12) 

o 

y(8)  ■ a(s)  + 0(s)  (13) 

*"Volume  cross  section"  is  a cross  section  per  unit  volume  of  scattering  or 
absorbing  centers  and  has  the  unit  of  inverse  length.  For  example,  a(s)  =* 
np(s)  aa( s)  where  n is  particle  number  density  and  o&  is  the  Mie  absorption 
cross  section.  In  Ref.  6,  these  quantities  were  referred  to  is 
"coefficients. " 


where  y(s)  Is  the  volume  extinction  cross  section. 


Equations  (1)  and  (2)  are  used  to  compute  the  combined  gas-plus-particle 
E/A  profiles  at  A^.  The  particle-only  equations  that  are  appropriate  for  the 
off-band  spectral  position  are  obtained  by  setting  k m o ( y = t s 1).  The 
results  for  particle-only  radiance  and  transmittance  are 

L 

np  " ^ V8)  V8)  JV8>  + Qs(8)}  ds  04) 


Q_  - J (s) 
T P 


J (s)  - a(8)  B (s) 
P p 


QS(s)  “ / JS(s,ft)  + 0T(s.n)]dft 


Jg(8,0)  ■ p(s,ft) 


QE(s,fi)  - Ta  [L(s,n)J  NE(sfO) 


L(s,n) 

QT(s,n)  - / T a(a)  Qt(o,s)  da 

o 


QT(a,s)  - J (a) 
1 P 


J Co)  * a(o)  B (a) 
P P 


-In  t ■ / y(s)  ds. 
o 


The  gas-only  equations  are  obtained  by  setting 

(t  a t=  1)  . The  gas-only  radiance  and  transmittance  are 
<*  p 


N - / t(s)  (^(s)  as 


Qt(s)  - y(s)  J (s) 


J_(s)  - x(s)  B (8) 
g g 


-In  x m I k(s)  y(8)  ds. 
8 o 


2.2  Transformation  to  Cylindrical  Coordinates 


Thus  far,  the  radiation  formulation  has  been  developed  In  terms  of  the 
observation  llne-of-slght  variable  s.  For  the  cylindrical  geometry  of  an 
axisymmetrlc  plume,  the  radial  coordinate  r is  a more  natural  coordinate.  The 
mathematical  detail  of  carrying  out  the  transformation  from  s to  r is  reported 
in  Ref.  6,  and  the  geometry  for  the  transformation  is  shown  in  Fig.  3.  In 
effect,  equations  of  the  form 


F * / J(s)  G(s)  ds, 


where  the  integration  is  over  the  full  extent  of  a chord  located  a transverse 
distance  z from  the  cylinder  diameter,  and  J(s)  is  a function  only  of  radius, 
are  transformed  into 


rar 


(29) 


F(*)  - 2 / J(r)  G(z,r)  — j j-jt 

z (r  ~ z J L 

where  R Is  the  cylinder  radius,  and 

G(z,r)  --j  [G+(z,r)  + G~(z,r)].  (30) 

G+(z,r)  and  G~(z,r)  are  the  values  of  G at  the  two  intersections  (Fig.  3)  of 
the  chord  8 at  z with  the  circle  of  radius  r.  The  reason  for  the  transfor- 
mation to  the  form  of  Eq.  (29)  is  not  that  it  facilitates  the  calculation  of 
E/A  profiles  (computer  coding  for  the  coordinate  s is,  in  fact,  somewhat 
easier  and  more  straightforward  than  for  r,  and  was  employed  in  the  previous 
code  EAPROF) , but  that  this  form  is  required  by  the  inversion  procedure  u;  ed 
in  this  work  (Section  3)  to  obtain  the  source  function  J(r)  from  a given 
F(z).  For  predicting  E/A  profiles  from  known  plume  properties,  the  kernel 
G(z,r)  is  a known  function.  In  inversion,  it  is  iteratively  computed  from  the 
retrieved  plume  properties  of  the  previous  iteration. 

The  transformations  of  the  observation  line-of-sight  equations  [(l.e., 
Eqs.  (1),  (12),  (14),  (23),  (24),  and  (27)]  to  the  form  of  Eq.  (29)  are  sum- 
marized. in  Table  1.  The  factoring  of  the  integrands  of  these  equations  into 
the  product  from  J(r)  G(z,r)  was  based  on  the  fact  that  only  the 
variables  a,  3,  p(fl),  ic,  Bg,  and  Bp  (and  thus  also  Jg  and  Jp)  are  functions 
of  r alone,  and  was  performed  to  isolate  the  particle  source  functions  in  the 
particle-only  equations  and  the  gas  source  functions  in  the  gas-only  and  the 
coupled  gas-plus-particle  equations.  The  factoring  for  the  particle-only  and 
coupled  gas-plus-particle  equations  also  makes  the  kernel  function  an  explic  t 
function  of  the  source  function  J.  In  E/A  profile  prediction,  this  explicit 
dependence  of  G on  J is  of  inconsequential  concern,  and  in  inversion  it  Is 
handled  by  using  the  previous  iteration  result  for  J in  calculation  G. 


The  actual  procedures  for  computing  the  various  components  of  G have  been 

discussed  previously:  ir*(z,r)  and  y*(z,r)  In  Ref.  1 end  Qg*(z,r)  In  Ref.  6. 

The  particle  transmlttances  r (z,r)  and  t.  (z,r)  are  computed  in  the  same 

a p 

manner  as  the  gas  transmlttances  t (z,r). 

2.3  Laser  Scattering  Equation 

The  scattered  radiance  at  angle  6 from  a beam  traversing  the  chord  s for 
the  case  where  only  single-scattering  is  important  is* 

I . L s L(s,0) 

1(e)  * ~2  I fl(s)  P(8,e)  exp  [-/  y(s')  ds'  - / y(o)  da]  ds  (37) 

where  IQ  is  the  Incident  beam  flux,  A is  the  sensor  entrance  aperture  area,  0 
is  the  distance  from  the  plume  axis  to  the  sensor  (it  is  assumed  that  D > >R), 
p is  the  scattering  phase  function,  and  o is  the  scattering  path  coordinate. 
The  intergration  over  s indicates  that  the  detector  field  of  view  encomposes 
the  entire  plume  diameter.  The  exponential  factor  accounts  for  extinction 
along  the  path  from  s « o to  s ■ 8 (o  » o)  and  then  on  to  a ■ L(s,0).  With  the 
definitions 

j(s,0)  -illl  P(6,0),  (38) 

s L(s,9) 

g(s , 0)  - exp  {-/  y(s')  ds'  - / Y (o)  da)]  (39) 

o o 


The  geometry  of  Fig.  3 is  applicable  to  the  laser  scattering  problem 
by  interpreting  s as  the  incident  axis  for  radiation,  O as  the  exit 
scattering  line  of  sight,  and  considering  scattering  only  within  the 
plane  of  the  figure  (0*  7T/ 2). 


Eq.  (37)  can  be  written  as 


J 


L 

f(0)  “ / J(s,e)  g(8,e)  ds.  (41) 

o 

As  for  the  E/A  equations,  this  result  is  now  transformed  to  cylindrical  coor- 
dinates. The  result  is 


f (z,0)  ■ 2 


R 

/ J(r,0)  G(z,r,6) 
z 


rdr 


2 \lk 


- *) 


(42) 


where 


J(r,6)  p(r,6),  (43) 

G(z,r,8)  - [g+(z,r,0)  + g (z,r,0)j,  (44) 

and  g*(z,r,0)  are  the  transmission  functions  that  correspond  to  the  two  scat- 
tering paths  from  i ■ o to  the  far  and  near  intersections  of  s with  r and  then 
on  to  the  boundary  at  angle  6.  The  method  of  calculation  of  g*(z,r,0)  is  a 
straightforward  numerical  integration  along  the  two  legs  of  the  scattering 
path  with  interpolation  for  y at  each  point  from  the  radial  profile  y(r). 


■ » 

r ■ 


i 
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3.  RETRIEVAL  PROCEDURES 


3.1  Introduction 

In  this  section,  the  methods  for  retrieving  the  plume  flow-field  proper- 
ties from  measured  (or  simulated)  E/A  and  laser  scattering  data  are  con- 
sidered. The  basic  method  of  retrieval  is  the  application  of  the  simple  Abel 
and  the  interative  Abel  inversion  procedures  derived  in  Ref.  1 to  a specific 
transverse  data  profile  to  retrieve  a specific  radial  profile.  The  procedure 
is  described  in  detail  in  Ref.  1 and  outlined  here  for  convenience. 

In  general,  we  have  an  equation  of  the  form 

R . 

F(z)  - 2 / J(r)  G{z,r,  J(r)J  (*5) 

where  F(z)  is  a known  transverse  profile  of  some  quantity  and  J(r)  is  the 
desired  result.  The  kernel  function  G is  a known,  calculable  function.  Its 
dependence  on  the  unknown  function  J(r)  is  handled  by  iteration  in  that  G is 
calculated  from  the  J(r)  of  the  preceding  iteration.  The  method  of  solution 
Is  outlined  in  Fig.  4 as  a calculation  flow  diagram.  The  function  F(z)  is 
entered  and  an  immediate  Abel  inversion  is  made  to  get  a first  guess  for 
J(r).  This  first  guess  would  be  the  desired  solution  if  G = 1 and  is 

<46> 


where 


H(r) 


2 / F(z) 


zdz 


z - r 


(47) 
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A ■ ABEL  INVERSION 


ENTER  F(z) 


INITIATING 

INVERSION 


ITERATION 

INVERSION 


CONVERGENCE 

TESTING 


Iterative  Abel  Inversion  Algorithm. 


With  this  first  guess,  * new  transverse  profile  ?'(*)  is  calculated  and  sub- 
tracted froa  the  true  profile.  An  Abel  Inversion  is  then  Bade  on  the  differ- 
ence profile  to  obtain  a correction  tern  to  J(r).  If  the  correction  tern  Is 
less  than  sone  convergence  criterion  e,  the  inversion  is  complete.  Otherwise, 
the  correction  is  made  to  J(r)  and  the  process  is  repeated  by  calculating  a 
new  transverse  profile  F'Ct)  with  the  new  J(r). 

The  input  transverse  profiles  required  by  the  inversion  procedure  are  ( is 

described  in  Section  1.3)  as  follows: 

il(z)  coupled  gas/particle  radiance  at 

t (a)  coupled  gae/particle  transmittance  at  Aj 

N (s)  partlcle-onlv  radiance  at  A, 

P ' 2 

r^(z)  particie-only  transmittance  at  A, 
t'(r,S)  scattering  efficiency  at  A^ 

T>«e  first  step  in  analysis  is  an  Abel  inversion  of  the  particle-only  transmit- 
tance profile  t (*)  to  get  the  radial  profile  of  the  volume  extinction  cro^s 
P 

section  r).  Next,  an  iterative  Abel  inversion  ia  made  of  f (2,8)  for  ea>  h 
(i.  The  radial  source  function  obtained  from  the  inversion  s 
•Hr, 3)  « 8(r)p(r,0)/A'v  where  0 < r > le  the  volume  scattering  cross  section, 
and  p(r,0)  is  the  scattering  phase  tuner f no.  An  integration  of  J(r,6),  over 
4?  steradians  for  each  r gives  U(r).  The  phase  function  is  then  obtain- d 
from  p(r,A)  * 4*  J(r ,li)/?(r).  and  the  volume  absorption  cross  section  -y 
*{r)  » v(r)  - Bit). 

At  this  point,  the  optical  properties  nt  the  particles  (except  for  their 
temperature)  are  completely  determined  in  so  far  as  they  are  needed  for  radia- 
tion calculations.  Note  also  that  nothing  has  been  assumed  about  the  parti- 
cles except  that  if  they  are  Irregularly  shaped,  they  are  randomly  oriented* 
If  assumptions  are  now  made  shout  what  the  particles  are  «wde  of,  what  their 


index  of  refraction  is,  and  that  they  are  spherical,  the  results  [that  is 

o(r),  0(r)  and  p(r,0)]  could  be  further  analyzed  to  yield  the  spatial  and 

1 

size  distributions  of  the  particles.  This  problem  is  not  considered  here. 

/ 

TJhe  next  step  in  the  analysis  is  an  iterative  Abel  inversion  of  the 

✓ 

particle-only  radiance  Np(z)  for  the  particle  temperature  Tp(r).  Now,  all  of 
the  particle  characteristics  are  determined,  and  the  final  step  can  be  per- 
formed (with  iterative  Abel  inversion)  by  solving  the  coupled  gas/particle  E/A 
profiles  N(z)  and  t(z)  for  the  gas  temperature  and  concentration  profiles. 

The  order  in  which  these  steps  are  performed  is  important  since  each  step 
requires  the  results  obtained  in  the  preceding  steps.  The  end  results  of  the 
inversions  are: 

T (r)  Gas  temperature 

c (r)  Gas  concentration 

g 

Tp(r)  Particle  temperature 

a(r)  Particle  volume  absorption  cross  section 

f5(r)  Particle  volume  scattering  cross  section 

p(r,0)  Particle  scattering  phase  function 

The  following  sections  give  the  steps  of  the  analyses  in  greater  detail. 

3.2  Retrieval  of  y(r) 

Of  all  the  equations  of  Table  1 that  define  the  transverse  profiles  in 
terms  of  radial  profiles,  Eq.  (34),  for  the  particle-only  transmittance,  is 
unique  because  the  kernel  function  G is  identically  unity.  Thus,  an  immediate 
Abel  inversion  of  -In  t (z)  yields  the  radial  volume  extinction  cross  section 


Transfer  Equations  in  Cylindrical  Coordinates 


8) 


The  second  step  in  the  retrieval  Involves  an  iterative  Abel  inversion 
of  the  laser  scattering  function  f(z,0)  for  each  0 to  obtain  J(r,0)  =» 
0(r)  p(r,0)/4ir  [see  Eqs.  (42)  and  (43)].  That  is 

J(r,0)  - . u [f(z,0)J  (49) 

where  IA  is  the  iterative  Abel  inversion  operator  defined  by  the  procedure  of 
Fig.  4.  The  iteration  initiation  assumption  G ■ 1 is  equivalent  to  the  as- 
sumption that  the  plume  is  nonscattering.  Integration  of  this  result  over  all 
solid  angles  yields 

/ J(r,0)  dft  * / p(r,0)  dfi  - 0(r)  (50) 

4n  **  4tt 

since  the  phase  function  is  normalized  to  4tt.  With  the  volume  scattering 
cross  section  0(r)  now  known,  the  volume  attenuation  cross  section  a(r)  can  be 
obtained  from 

o(r)  - y(r)  - 0(r),  (5!) 

and  the  scattering  phase  function  from 


For  some  applications,  it  nay  be  reasonable  to  assuae  that 

a,  8,  and  p(0)  to  not  depend  on  radius.  In  this  case,  the  line-of-sight 

equation  (Eq.  (42)]  reduces  to 

R 

f(z,0)  -j-p(o)  2 f G(z,r, 6)  — 

z (r  - z J 72 

and  $p(0)  can  be  obtained  directly  from  data  taken  for  a single  transverse 

value  z (presumably  z * o for  best  signal)  without  inversion  by 

J(0)  -|^.p(e) g . (54) 

2 / G(z,r,9) 

z [r  - z J z 

As  for  the  general  case,  an  integration  over  all  solid  angles  allows  for  a 
determination  of  8 (and  thus  a)  and  p(0)  individually.  The  consequences  of 
this  simplification  to  the  diagnostic  itself  are  not  great  since  the  iterative 
Abel  inversion  for  the  general  case  Is  easy  to  implement  and  is  fast.  A 
significant  simplification  to  experimental  measurement  would  be  effected, 
however,  since  now  the  laser  source  would  not  have  to  be  traversed  in  z. 

3.4  Retrieval  of  Tp(r) 

The  particle  temperature  profile  is  retrieved  by  iterative  Abel  inversion 
on  the  particle-only  radiance  profile.  The  retrieve/  source  function  Is  (see 
Eq.  (33),  Table  lj. 


a(r)  B (r)  - IA  [N  (z)J. 
P P 


(55) 
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With  a(r)  known,  B,  rj  and  thus  Tp(r)  can  be  found.  The  iteration  initiation 
assumption  G * 1 for  this  inversion  is  equivalent  to  the  assumption  that  the 
plume  is  nonscattering. 

At  this  point,  all  of  the  particle  properties  that  are  required  for 
further  analysis,  or  for  radiation  calculations  in  general,  have  been  ob- 
tained. 


3.5  Retrieved  of  Tg(r)  and  ic(r) 

The  retrieval  for  gas  temperature  Tg  and  volume  gas  absorption  cross 
section  k proceeds  by  a simultaneous  iterative  Abel  inversion  on  the  total 
gas-particle  emission  and  absorption  profiles  [Eqs.  (31)  and  (32),  Table  lj. 
The  need  for  a simultaneous  solution  is  that  the  emission  and  transmission 
profiles  are  coupled  through  the  dependence  of  k on  T , This  coupling  did  not 

O 

occur  in  the  particle-only  retrieval  of  Tp  and  y because  y can  be  obtained  in 
a single  Abel  inversion.  The  retrieved  source  functions  are 


<(r)  B (r) 

o 

ic(r) 


- IA 


N(z) 

x(z) 


(56) 


The  retrieval  of  the  gas  temperature  Tg(r)  and  the  gas  concentration  cg(r) 
from  these  two  source  functions  is  described  in  Ref.  1.  The  initiation  as- 
sumption G - 1 for  the  iterative  Abel  inversion  is  equivalent  to  the  assump- 
tion that  the  source  is  a thin,  gas-only  plume. 
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4.  APPLICATION 


The  particle**  and  gas-property  retrieval  schemes  developed  In  Section  3 
were  applied  to  the  minimum  smoke  propellant  (MSP),  advanced  liquid  propellant 
(ALP),  and  reduced  smoke  propellant  (RSP)  models  considered  in  the  previous 
phase  of  this  work  (Ref.  6).  In  Section  4.1,  the  flow-field  properties  of  the 
three  plume  models  are  summarized,  and  in  Section  4.2,  the  computed  transverse 
E/A  and  laser  scattering  profiles  are  discussed.  In  Sections  4.3  and  4.4, 
respectively,  the  particle  and  gas  retrieval  results  are  presented. 

4. 1 Plume  Model  Summary 

A summary  of  the  plume  model  parameters  is  given  in  Table  2 and  the 
radial  temperature  and  concentrations  are  shown  in  Figs.  5-7.  The  differen- 
tial scattering  cross  sections  used  for  the  three  particle  species  A^O-j,  C, 
and  Zr02  are  shown  in  Figs.  8-10.  These  scattering  cross  sections  (and  the 
cross  sections  in  Table  2)  for  AI2O3  and  Zr0£  were  obtained  for  the  bimodal 
particle  size  distribution  described  in  Ref.  6.  The  peaks  of  the  distribution 
are  at  particle  radii  of  0.05  and  0.5  pm.  The  index  of  refraction  used  was 

m “ 2.0  - 0.01  i at  X * 2.51  pm  for  AI2O3  and  m ■ 1.2  - 0.01  i at  X * 3.38  pm 

for  Zr02#  The  8cattering  cross  sections  for  C were  obtained  with  a size 
distribution  function  that  peaks  at  a radius  of  0.05  pm  and  falls  by  two 

orders  of  magnitude  by  0.5  pm.  The  index  of  refraction  used  was  m ■ 2.6  - 1.0 

i.  The  gas  band  model  parameters  for  B^O  at  2.51  pm  and  HC1  at  3.89  pm  are 

discussed  in  Ref.  6. 
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Table  2.  Summary  of  Plume  Model  Data, 


Model 

MSP 

ALP 

RSP(a) (b) 

X(pm) 

2.51 

3.89 

2.51 

Gas 

h2o 

Ha 

h2o 

Particle 

*l2o3<W 

C 

Zr02 

R(cm) 

10 

7 

5.69 

c (mole  fraction) 

0.15 

0.117 

0.155 

p (atm) 

1.0 

1.18 

— 

—3 

np  (cm  J) 

105 

1.37  x 108 

— 

„ t 2' 

a (cm  j 
a 

3.20  x 10'9 

4.62  x 10"10 

1.41  x 10"9 

a (cm^) 
s 

5.86  x 10"8 

2.17  x 10"10 

8.41  x 10"9 

Lj  (cm) 

3.0 

2.5 

3.0 

L2  (cm) 

57 

83 

44 

T„<K) 

1460 

eN 

0.75 

0.25 

0.25 

X wavelength 

R plume  radius 

c gas  concentration 

p total  gas  pressure 

n particle  loading 

0^  absorption  cross  section 

a scattering  cross  section 

Lj  exit  plane  - observation  plane  distance 

L<2  observation  plane  - end  plane  distance 
nozzle  temperature 
eN  nozzle  emlaslvity 

(a)  Low-temperature  model  of  Ref.  6. 

(b)  With  Imaginary  Index  k *=  0.01. 


4.2  Emission/ Attenuation  Profiles  and  Laser  Scattering  Results 

For  all  three  cases,  transverse  profile  calculations  were  performed  for  a 
Lorentz  lineshape  and  the  Curtis-Godson  nonuniformity  approximation.  For 
numerical  integration  along  the  observation  lines  of  sight,  the  plume  was 
divided  into  10  annular  zones.  Integrations  along  scattering  lines  of  sight 
employed  a grid  of  10  equal-size  intervals.  Angular  integration  was  performed 
with  the  11-point  scattering  angle  grid  0 * 0,  5,  15,  25,  35,  45,  60,  90,  120, 
150  and  180°  and  an  azimuthal  grid  of  16  equal-size  intervals  (A<j>  * 22.5°). 
The  laser  scattering  function  f(z,0)  was  computed  on  the  8-point  grid  0 ° 0, 

Li 

10,  20,  30,  45,  90,  135,  and  180°. 

The  computed  transverse  profiles  for  total  gas-particle  radiance  and 
transmittance  are  shown  in  Figs.  11-16,  and  duplicate  the  results  of  Ref.  6. 
The  MSP  and  RSP  models  represent  plumes  whose  radiation  and  attenuation  are, 
for  the  most  part,  dominated  by  the  gas  component  of  the  plume.  In  the  ALP 
model,  however,  particles  play  the  dominant  role.  The  loading  of  the  ALP 
moiel  is  so  large,  in  fact,  that  the  condition  of  signal-scattering  required 
by  the  radiation  model  is  most  assuredly  violated,  and  neither  the  radiation 
nor  the  extinction  results  for  this  model  should  be  believed.  However,  the 
results  can  be  used  to  test  the  inversion  procedures  because  the  same  single- 
scattering model  is  employed  in  the  inversion. 

Representative  scattering  efficiency  function  results  are  given  in  Figs. 
17-22.  The  angular  variation  of  these  functions  essentially  reflect  the 
angular  variation  of  the  particle  scattering  phase  functions  (Figs.  8-10). 
(Note,  however,  that  these  functions  are  not  required  to  be  symmetric  at 
0 n 0 and  180°,  while  the  phase  functions  are.)  The  overall  decrease  in 
magnitude  of  the  angular  variation  with  increasing  z is  caused  by  the  de- 
creased scattering  path  length  and  the  increased  amount  of  plume  that  the 
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scattered  radiation  wist  pass  through  before  reaching  the  sensor.  In  the  RSP 
model  (Fig.  21),  the  decrease  - ■ - ^nitude  is  further  enhanced  by  the  rapid 

decrease  in  particle  density  with  radius  (Fig.  7). 

The  E/A  and  scattering  results  of  Figs.  11-22  are  the  simulated  experi- 
mental data  to  which  the  inversion  algorithms  are  applied  in  the  following  two 
sections. 

4.3  Particle  Retrieval  Results 

For  each  of  the  three  models,  retrieval  was  made  for 

y(r),  8(r),  a(r),  p(r,0)  and  Tp(r)  as  described  in  Section  3.*  All  of  the 
computation  conditions  used  to  generate  the  transverse  profiles  (e.g. , angle 
grids  and  number  of  zones)  were  duplicated  in  the  inversions.  For  the  re- 
trieval of  p(r,0)  and  Tp(r),  which  require  iteration,  the  convergence  re- 
quirement imposed  was  that  the  radial  root-aean-square  difference  between 
successive  iterations  for  the  respective  source  functions  be  less  than  or 
equal  to  0.1%. 

In  all  cases  but  one,  the  retrieved  radial  profiles  and  phase  functions 
were  exact  to  within  the  convergence  criterion.  The  one  exception  was  the  RSP 
plume  model.  Here,  the  retrieval  was  exact  near  the  plume  center,  but  became 
increasingly  noisy  in  the  outer  zones  of  the  plume.  In  the  outermost  zone,  in 
fact,  retrieved  volume  cross  sections  became  negative,  and  retrieval  for  Tp 
was  idefinite.  The  retrieval  errors  are  shown  in  Fig.  23.  The  cause  ot  the 
errors  is  the  extreme  sparseness  of  particles  in  the  outer  regions  of  the 
plume  (Fig.  7),  and  as  such  is  not  particularly  significant.  The  large  errors 


For  comparison  with  input  date,  y*  ti  (<r  + ff_),  fim  n , Ct*  nD  <7a»  and 
p(0)  “ (4 tf/<rg)  d<7fl/d Q. 
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Pig.  23.  Particle  Properties  Retrieval  Errors 
for  the  RSP  Plune  Model 


i 


i 

| 

i 


L 


55 


merely  reflect  the  fact  that  we  cannot  recover  the  properties  of  nonexistent 
(or  nearly  so)  particles. 

4.4  Gas  Retrieval  Results 

For  each  of  the  three  plume  models,  retrieval  was  made  for  gas  tempera- 
ture and  concentration  using  the  gas-only  and  first-order,  off-band  correct  on 
procedures  described  in  Ref.  6,  and  the  fully-coupled  gas-plus-particle  pro- 
cedure developed  here  in  Section  3.  In  the  gas-only  procedure,  the  transverse 
profiles  of  N and  t were  treated  as  if  they  arose  from  a purely  gaseous  plume 
and  used  directly  as  input  in  a gas-only  inversion.  In  the  first-order,  off- 
band  correction  procedure,  a gas-only  inversion  was  made  on  the 

profiles  N - N and  t/t  . The  results  of  these  inversions  are  shown  in  Figs. 
P P 

24-29.  The  results  tabled  GAS  are  the  gas-only  inversion  results,  and  the 
results  labeled  FO/OB  are  those  obtained  with  the  first-order,  off-band  cor- 
rection procedure.  For  the  RSP  model,  the  FO/OB  results  are  indistinguishable 
from  the  true  profiles.  The  inversions  were  performed  with  the  same  computa- 
tion conditions  (e.g.  llneshape  and  number  of  zones)  used  to  generate  the 
transverse  profiles  and  with  the  convergence  criteria  that  the  radial  root- 
mean-square  difference  between  iteration  results  for  T and  c„  should  be  less 

o o 

than  or  equal  to  0.1Z.  These  inversion  results  duplicate  the  results  obtained 
in  Ref.  6. 

The  results  of  the  fully-coupled  inversion  are  not  shown  explicitly  in 
Figs.  24-29.  In  all  cases,  the  recovered  radial  profiles  were  within  the 
convergence  criteria  and  are  indistinguishable  from  the  true  profiles  when 
plotted. 
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C (mole  fraction) 


Fig.  29.  Gas  Concentration  Retrieval  Results  for 
the  RSP  Plume  Model 
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5.  SUMMARY  AND  CONCLUSIONS 


In  the  first  phase  of  this  work  (Ref.  6).  and  repeated  here,  it  was  shown 
that  even  for  plumes  as  lightly  loaded  with  particles  as  the  MSP  and  RSP 
models,  the  effects  of  particles  could  not  simply  be  ignored  in  attempting  to 
retrieve  gas  temperature  and  concentration  from  two-phase  plumes.  Retrieved 
temperatures  were  invariably  too  low,  and  concentrations  were  too  high.  With 
partial  account  of  particle  effects  through  the  use  of  the  first-order,  off- 
band  correction  procedure,  retrieval  for  these  two  lightly  loaded  plumes  was 
substantially  improved.  For  all  practical  purposes,  the  first-order  correc- 
tion procedure  suffices  as  an  accurate  diagnostic  for  the  gas  properties  of 
these  plumes.  For  the  more  heavily  loaded  ALP  plume  model,  however,  even  this 
correction  procedure  could  not  yield  good  gas  retrieval  results.  In  response 
to  this  problem,  the  fully-coupled  gas-plus-particle  retrieval  procedure 
described  in  this  report  was  developed.  This  new  procedure  is  capable  of 
arbitrarily  accurate  gas  property  retrieval  from  ideal  data.* 

In  addition  to  the  improved  gas  property  retrieval  procedure,  a meth- 
odology for  retrieving  particle  properties  from  appropriate  E/A  and  laser 
scattering  data  was  developed.  These  rocedures  were  also  shown  able  to 
retrieve  particle  properties  to  arbitrary  accuracy  from  Ideal  data. 

In  all,  the  results  of  this  work  demonstrate  that  a complete  diagnostics 
procedure  for  retrieving  both  particle  and  gas  properties  from  two-phase,  low- 
visibility  propellant  plumes  could  be  developed.  The  fundamentals  of  the 
procedure  have  been  developed  and  tested  here.  Additional  work  must  be  made 


A minor  penalty  of  substantially  increased  computer  running  time  is  paid 
by  using  the  fully  coupled  retrieval  procedure  over  either  the  gas-only 
or  first-order,  off-band  procedures. 
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on  error  propagation  effects,  and  procedures  for  retrieving  particle  size  and 
spatial  distributions  would  enhance  the  value  of  the  diagnostic. 
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